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cn . 1. Introduction 

r^ • Ever since Newton's formulation of particle mechanics over three hundred years ago, 

^N , one of the most widely used methods of modeling physical systems is via an initial value 

Jh I formulation. For particle mechanics, the idea is that one specifies the initial position yo 

^ • and initial momentum pQ of the particle, and one then determines the particle's path y{t) 

^ [ in time by solving the initial value problem consisting of Newton's equations (an ordinary 

C""^ ' differential equation system) 

^ im = p(m, 

\ together with the initial conditions y(*o) = Vq and p{tQ) = pQ. Analogously, to study the 

^7j ' vibrational motion ■ip{x,t) of a stretched string, one specifies the initial displacement il^o{x) 

(-—\ , and initial velocity I'd (2;) of the string, and one then determines the subsequent motion 

^O ' ''Pix^t) by solving the initial value problem consisting of the (first order) wave equation (a 

0^ . partial differential equation system) 

T^' ; dtip{x,t) = u{x,t), 

0| dtv{x,t) = adxxi^{x,t), 

^"^ \ together with the initial conditions ip{x,tQ) = 'ipo{x) and iy{x,tQ) = i'o(2^)- 

^ ' Although one of the signature properties of Einstein's theory of gravity is its spacetime 

covariant character, it too admits an initial value formulation. We recall that to model the 
gravitational interactions of physical systems, general relativity incorporates the physics 
C^_' of the gravitational field into the geometry of spacetime, which is specified by a four- 

dimensional manifold M^ together with a metric g of signature (—,+,+,+). Einstein's 
equations relate the spacetime metric to the non-gravitational "matter fields" ^ which are 
present in the spacetime by requiring that the Einstein tensor G^,^ := R^,^ — {l/2)Rg^y of 
g (here R^y is the Ricci tensor of g and R is its scalar curvature) be proportional to the 
stress-energy tensor T^y of ^ and g; specifically. 



X 



G^M = ^^T^u[gM, (1) 
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where k = ^ttQn with Qj^j Newton's gravitational constant. One of the most effective ways 
to construct and study a solution {M'^^g,^!) to Einstein's equations ([TJ is by specifying 
a set of initial data corresponding to an initial "state of the universe", and then using 
Einstein's equations (reformulated as an initial value problem) to evolve this data into the 
corresponding spacetime solution. This chapter discusses the details of this initial value 
formulation for Einstein's theory: the specific nature of the initial data sets for which the 
initial value problem works, how one finds such initial data sets, how one reformulates 
Einstein's equations ([1]) as an initial value problem and verifies that it is well-posed, and 
some open problems related to the evolution of initial data sets. 

One very important feature of the initial value formulation of Einstein's theory is that 
the initial data must satisfy a set of constraint equations. This is a familiar feature of 
Maxwell's theory of electromagnetism as well: An initial data set for Maxwell's theory 
consists of a pair of spatial vector fields {Eq (x) , Bq {x) ) which are required to satisfy the 
Maxwell constraint equations 

V • ^0 = 0, 

V • ilo = 

(presuming that the charge density is zero). A solution {E{x,t),B{x,t)) is obtained from 
the initial data via the evolution equationqj 

dtB{x,t) = V X E{x,t), 

dtE{x,t) = -V X B{x,t). 

It is important to note that presuming a data set {Eq{x),Bq{x)) satisfies the Maxwell 
constraints, it follows from the evolution equations that the solution generated from this 
data set satisfies the constraints V • E{x,t) = and V • B{x,t) = for all times t. 

To understand the initial value formulation of Einstein's theory and how it produces 
spacetimes (M^, (7, ^) which satisfy Einstein's equations, it is useful to start by examining 
spatial foliations of such a spacetime, with a focus on the geometric fields induced on the 
spatial leaves of such a foliation and how they relate to the spacetime metric g. We do 
this in Section [2j Included in this section is a discussion of the Gauss-Codazzi-Mainardi 
equations which relate the spacetime curvature of g to the induced fields. Based on these 
equations, we derive the constraint and evolution equations which comprise the initial value 
formulation of Einstein's theory, and then in Section [3] we discuss the proof that this initial 
value problem is well-posed. Next, we examine methods for obtaining initial data sets 
which satisfy the Einstein constraint equations. In Section [5] we focus on the conformal 
method, in Section [5] we discuss the closely related conformal thin sandwich method, and in 
Section [6l we explore how such solutions can be obtained via gluing. Our discussion of the 
Einstein evolution equations is very brief; in Section [7] we comment primarily on globally 
hyperbolic solutions with singularities and the strong cosmic censorship conjecture, along 
with brief remarks concerning the stability of Minkowski spacetime, the Kerr solutions, 
and certain expanding cosmological solutions of the Einstein equations. 



Here and throughout this chapter, we set c, the speed of hght equal to one. 
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2. (3 + l)-FOLIATIONS AND THE DERIVATION OF THE EiNSTEIN CONSTRAINT AND 

EVOLUTION EQUATIONS 

To see how to build spacetime solutions of the Einstein field equations from initial 
data on a Riemannian manifold, it is useful to start with such a spacetime (M^,(7, ^), 
specify a spatial foliation of that spacetime (presuming that it admits suclo), determine 
the geometric quantities which are well-defined on the leaves of that foliation, and then 
calculate what the Einstein equations ^ tell us about these geometric quantities. We do 
this here. 

A three-dimensional hypersurface S^, smoothl}o embedded in {M'^,g) (with embedding 
map z : S^ — 7- M^), is defined to be spacelike (or spatial) if the induced bilinear form 
7 := i*g on S^ is Riemannian. Equivalently, the embedded hypersurface is spacelike if all 
vectors V tangent to i(S^) at p are spacelike with respect to g (so that g{V,V) > 0). A 
spatial foliation of the spacetime {M^,g) is a smooth one-parameter family ij : S^ — )• M^ 
(for t G M) of embedded spacelike hypersurfaces such that each point p £ M^ is contained 
in one and only one of the hypersurfaces ii(S^) (referred to as the "leaves" of the foliation). 
We presume that it is also true that for each p G M^, there is one and only one point q E S'^ 
such that ifiq) = p for some t. 

We focus now on a particular leaf «i(S^) of a chosen foliation it of a chosen spacetime 
{M'^jg). The definition of spatial foliation (above) guarantees that ii(S'^) is equipped 
with a Riemannian metric ^m, which consequently defines a covariant derivative Vm and 
the corresponding Riemann, Ricci, and scalar curvature quantities, all intrinsic to ft(S^). 
In addition to specifying these intrinsic geometric quantities, the foliation if equips it(S^) 
with extrinsic geometric quantities, including (i) e_|_w , a future-pointing unit-length timelike 
vector field, orthogonal to all vectors tangent to it{T,^); (ii) ^m, a unit-length one-form field 
such that (at each point in it{T,^)), {Ouu^±[t]) = 1 ^^'^ (^"""(p)?^) = for all vectors V 
tangent to it(S^); and Km, the second fundamental form, defined by 

K^t0,V) := g{Due^[t^,V), (2) 

where D denotes the spacetime covariant derivative and where U and V are vector fields 
tangent to ii(S^). It is straightforward to show that K^^-^ is a well-defined (spatial) tensor 
field with respect to the tangent spaces of ii(S^), and further that it is symmetric. By con- 
trast, e_|_w and ^^i are not spatially tensorial; rather, they act tensorially on the spacetime 

tangent spaces of M^, restricted to the subset ii(S'^) C M^. 

To relate the foliation-related quantities 7[i],V[t], the spatial curvatures, -f^[t], e_|_[(] and 
0m to the spacetime geometry, it is useful to specify foliation-compatible bases. We obtain a 
(local in space) foliation-compatible coordinate basis by choosing coordinates {xa}\(a=i,2,3) 
locally on S^, using the foliation to transport these to its leaves it(S^), and adding the 



So long as the spacetime is globally hyperbolic, as defined below, it admits a spatial foliation. Spacetimes 



which are not globally hyperbolic do not generally admit spatial foliations 

To simplify the discussion and avoid issues of regularity, we 
smooth, by which we mean that they are infinitely diflerentiable 



To simplify the discussion and avoid issues of regularity, we presume that all manifolds and fields are 
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parameter t to produce the spacetime coordinates {xa,t} =: {xa}\ (0=1,2,3,0) ■ The corre- 
sponding dual coordinate bases {da} and {dx°'} are fohation-compatible in the sense that 
dt annihilates vectors which are tangent to the leaves, dt is transverse to the leaves, and 
da are tangent to the leaves; however, we note that dt and dt need not be timelikqj ev- 
erywhere. Alternatively, one may choose the foliation-compatible basis {enfuda} and its 
dual {^M;^M := dx"' + M^dt}; here M^i are the components of the shift vector, which is 
the spatial projection of dt- 

dt = N[t]e^[t] + M^^^da. (3) 

The scalar Nm is called the lapse function. 

On the leaf it{T,'^), one easily verifies that the spacetime metric g is related to the spatial 
metric ju] by the identity (and its expansion) 

9 = 7[t]-^[[]^[[] (4) 

= l[t]ab{dx'' + Mli]dt){dx' + Ml^dt) - Nl^dt\ (5) 

The expression relating the spacetime covariant derivative D (compatible with g) and the 
spatial covariant derivative Vrji (compatible with 7r^i) is not so easily derived; however it 
does follow from the definition (j2]) of the second fundamental form that for any pair of 
vectors U and V tangent to ij(Il^), one has 

DuV = V^t]uy + K^t]{U,V)e^^t], (6) 

thus identifying Ku] with the surface- normal projection of the spacetime covariant deriv- 
ative. One also verifies from this definition, as well as from the Lie derivative identity 
C-YQaP = DaYp + DpYa, that 

^dtl[t]ab = '^^[t]K[t]ab + ^Myt]1[t\ab- (7) 

Notably, this result involves not only quantities on the slice iw(S^), but also the rate of 
change of one of these quantities-7w -along the foliation. 

It remains to consider the relationship between the spacetime curvature and the spa- 
tial curvatures. Since the calculations which lead to formulas for the components of the 
spacetime curvature tensor in terms of components of the spatial curvature tensor along 
with other quantities intrinsic to leaves of a chosen foliation are analogous to analyses first 
done (for Riemannian rather than for Lorentzian ambient geometries) by Gauss, Codazzi 
and later Mainardi, these formulas are often referred to as the Gauss-Codazzi-Mainardi 
equations. Using 

7^°^^5 := (r , D,^D,,ep - D,,D,^ep - D[e„e,]e;3> (8) 

to represent the spacetime curvature tensor (with respect to spacetime dual bases {ca} 
and {6^}) and analogously 

R\cd ■■= {^^ Ve,Ve,efe - Ve,Ve,eb - V [e.,e,]ef,) (9) 



dt is timelike if g{dt,dt) < and it is null if g{dt,dt) < 0; similarly dt is timelike or null or spacelike 
depending on the sign of g~^{dt,dt). 





'^\cd - R%cd + K^Km - 


- K^Khc, 






n\,^ = Ve^KS-Ve,Kc, 


T^\i>± -- 


= -Ce^K^ + K-"'K,r^b 


+ ^ V^V.iV, 
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to represent the spatial curvature tensor (with respect to spatial dual bases {cq} and {0^}), 
we find that these equations take the following fornu: 

(10) 

(11) 
and 

^ (12) 

Note that, to avoid notational clutter, in equations (f9])- (|12p we have removed the "[t]" 
subscripts from foliation-defined quantities such as V[t] and K^t^ab- 

If the spacetime {M^,g, ^) under consideration satisfies the Einstein field equations ([1]) 
and if a spatial foliation it : E'^ — )• M^ has been specified, what do these equations tell 
us about the foliation-defined geometric quantities 7m and Ku-t? Combining the definition 
of the Einstein tensor Ga/s = Tl-ap ~ 2'^9ai3 with the Gauss-Codazzi-Mainardi equations 
([9|)-(|12p. we readily obtain the following 

= ^Gx± = fr^±[7,V',^], (13) 

'c) = G^d = KTUl.^,A, (14) 

K',+K',Kab + ^VaybN = Gab-\9abGl (15) 

= KTabil, ^Jj, vr] - -gabT^b, ^/J, tt], 

where {ip, vr) represent the non-gravitational fields and their derivatives, after they have 
been decomposed into quantities which are well-defined with respect to the leaves of the 
foliationO 

The first two sets of these equations (fT3]) - (fH]) have the very interesting feature that 
(presuming the quantity T_i±[y , tp , tt] behaves well), they involve quantities (7[t] and K[^, 
not e^) defined as tensors on the leaf it(S^), and do not directly involve any time derivatives 
of these quantities. By contrast, the remaining set of these equations ([7]) does involve the 
time derivative of Km . 

There is nothing special about it(S^) or any other leaf of the foliation. Nor is there 
anything special about the (arbitrarily) chosen foliation. Hence we see that on any spacelike 
hypersurface embedded in a spacetime (M^,g', ^) which satisfies the Einstein equations, 
the surface quantities {'y , K , ip , it) must satisfy (fT3l) - (fT4l) . These are hence known as the 
initial value constraint equations for Einstein's theory. 

The initial value problem (or Cauchy problem) for Einstein's theory turns things around. 
It addresses the following question: If we fix a three-dimensional manifold S^ and if 



R + K.dK'"'' - 


[K? - 


^cK^i- 


- ^diK 


^e_i_Kab + Rab 


- 2Kac 



The most straightforward method of deriving the Gauss-Codazzi-Mainardi equations is by combining 
the expression (|6]) relating D, V and K with the expressions (|8]) and (|9} for the curvatures. Details of these 
calculations appear in [43) . 

For Maxwell's electromagnetic field, (i/',7r) correspond to {Ba,Ea)- 
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we choose a Riemannian metric 7^6 and a symmetric tensor K^d (and possibly, non- 
gravitational field data if) and vr on S^ as well) which satisfy the constraint equations 
(|13p - (|14p . is it always true that there exists a spacetime (M^,g(,^) which satisfies the 
(spacetime) Einstein field equations ([1]), and which contains an embedded submanifold 
z(S^) C M^ for which 7^^ is the induced metric i*gah and K^d is the induced second fun- 
damental form? Further, if such a solution exists for a given set of initial data, is it unique 
in any sense? We address this question in the next section. 

3. Well-posedness of the initial value problem for Einstein's equations 

It appears, from our discussion in Section [2l that the construction of a (vacuum) solution 
of the Einstein equations from initial data is a straightforward enterprise (especially with 
a computer readily at hand): a) On a specified three-dimensional manifold S^, we choose 
an initial data set (7ab,-f^cd) which satisfies the (vacuum) constraint equations 

R + K.dK^'' - {Klf = (16) 

V,K',-Vd{K) = ^; (17) 

b) we freely choose the lapse as a one-parameter (t) family of positive scalar functions 
N{x,t) on S^, and the shift as a one-parameter family of vector fields M^{x,t) on S^; c) 
using the evolution equations 

^dtlab = 2NK[t]ab + ^Anab, (18) 

J^a^Kab = -NRab + 2NKacK- NK'^Kab-Va^bN + CuKab, (19) 

we evolve the initial data set into a one-parameter family of Riemannian metrics jab{x,t) 
and symmetric tensors Kab{x,t); and d) we construct the Lorentz metric 

9 = lab{x,t){dx'' + M''{x,t)dt){dx^ + M\x,t)dt) - N^{x,t)dt^ (20) 

on the spacetime manifold Af^ = S^ x M and verify that it is a solution of the vacuum 
Einstein equations Ga/s = 0. 

To determine that this construction procedure in fact works, at least for some choices of 
the lapse and shift, one needs to prove a well-posedness theorem. Such a theorem was first 
proven for Einstein's theory during the 1950's by Yvonne Choquet-Bruhat [29]. This early 
work leaves the issue of uniquenness unsettled to a certain extent. However, the later work 
of Choquet-Bruhat with Robert Geroch |13) provides a strong form of uniqueness. Before 
stating this combined result, we find it useful to establish some terminology: 

A spacetime {M'^,g) is globally hyperbolic if there exists an embedded spatial hypersur- 
face i(S^) C M^ such that every future and past inextendible causal patljj intersects i(S'^) 
once and only once. If such a hypersurface exists, it is called a Cauchy surface. A spacetime 
is a globally hyperbolic development ( "gh-development" ) of a specified set of initial data 



A smooth path r) : I -^ A/* is causal if its tangent vector V is always either timehke {giV, V) < 0) or 
null (g(V,V) = 0). Such a path is inextendible if there does not exist a smooth path i) which contains rj as 
a proper subset. 
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(7, K) on S'^ if (i) it as a solution of the Einstein field equations; (ii) it is globally hyper- 
bolic; and (iii) there exists an embedded hypersurface ?(S^) such that {M^,g) induces the 
specified initial data (7,-?^) on «(S^) in the sense that 7 = i*g and an appropriately modi- 
fied version of ([2]) holds. A spacetime (M^,^) is a maximal globally hyperbolic developm,ent 
of (7, K) on S^ if every other gh-development of the same set of data is diffeomorphic to 
a subset of {M^,g). 

Using this terminology, we have the key result [29^ [T3] : 

Theorem 1 (Well-posedness of Einstein's vacuum equations). For any smooth set of initial 
data (7, -ftT) on Y? which satisfies the vacuum constraint equations (fT6]l - (fT7|l . there exists a 
unique (up to diffeomorphism) maximal globally hyperbolic development. 

Before discussing how to prove this result and discussing ways in which the result can 
be generalized, we note a number of its features: First, Theorem [1] is fundamentally a local 
existence and uniqueness result. It guarantees that for any initial data set, there exist 
spacetimes which are gh-developments of that data, and there exists a unique maximal 
gh-development, but it says nothing regarding whether that maximal development lasts 
for a finite or an infinite amount of proper (or coordinate) time. Second, Theorem [J 
guarantees that there are globally hyperbolic spacetimes evolving from any given set of 
initial data, but it says nothing about non-globally hyperbolic spacetimes containing an 
embedded hypersurface with that data. In many cases — the Taub-NUT spacetime j33] and 
the various generalized Taub-NUT spacetimes ^50] are prominent examples — the maximal 
development of a given set of data can be extended smoothly across a Cauchy horizon 
to a non-globally hyperbolic spacetime solution. A Cauchy horizon in a given spacetime 
{M'^,g) {not to be confused with a Cauchy surface, as discussed above) is a hypersurface 
71 embedded in (M^,gf) which is null (i.e., at each point ofTi, there is a vector tangent to 
Ti which is null) and which lies on the boundary between the region in [M'^^g) which is 
globally hyperbolic and the region which is not. The initial value problem tells us essentially 
nothing about whether or not such extensions exist for a given set of initial data. Third, 
Theorem [1] guarantees that there are choices of the lapse and shift which can be used in 
evolving initial data sets to produce gh-developments (numerical or otherwise), but it does 
not tell us whether this is true just for certain particular choices, or for all such choices. 
The proof of Theorem [H as we see below, does provide some information regarding this 
issue, but does not resolve it. Finally, we note that while well-posedness theorems usually 
contain statements regarding continuity of the map from the domain space of initial data 
sets to the range space of solutions, this is not true of Theorem [H as stated. Such a 
result does hold in an appropriate form for Einstein's equations; however, to avoid the 
unenlightening detail needed to state it, we do not include it here. We refer the interested 
reader to Chapter 15 of [55]. 

The key for proving that the initial value problem for Einstein's field equations is well- 
posed is to show that the system ([T| can be, in a certain sense, transformed into a hyperbolic 
PDE systeno for which well-posedness is well-established. The vacuum Einstein equation 



A PDE system is hyperbolic if, roughly speaking, it has the characteristics of a system of wave 
equations. A precise definition of hyperboUc PDE systems is found, for example, in [27) . and also in [55| . 
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system, which we can write as TZa/s = 0, is not itself hyperbohc. The system 

7^a/3 + -Czgal3 = 0, (21) 

with 

Za ■■= -g'^idug^a - Idagut.) + J'aig] (22) 

for J-'a[g] a specified function of ga/s but not of its derivatives, takes the form 

g''^d,d^g^^ = Qap[g,dg], (23) 

and therefore is a hyperbohc PDE system for the metric ga/s- This does not lead to the 
Einstein system itself being hyperbolic. However, using the Bianchi identities, one can 
show that the vector field Z"' satisfies the homogeneous hyperbolic system 

g'^'d.d^Z'' = -n$Z^, (24) 

so that the combined system of (|2i p -(|24 p is hyperbolic. Further, one can show (see, for 
example. Chapter 14 of [55], or Section 3.1 of [8j) that, so long as the geometric initial data 
(7, i^) satisfy the constraints, one can choose the corresponding initial values of gap and of 
dtgajB for the system (l2ip -(fM|l in such a way that the initial values of Z° and dtZ'^ vanish. 
It then follows from standard results regarding (nonlinear) hyperbolic systems (see, for 
example. Chapter 9 of [55]) that the system (f2T]l - (pl]l is well-posed; further, it follows that 
the solution {gap{x,t),Z°'{x,t)) has Z°'{x,t) vanishing for all time. Hence {gapix,t),0) is 
a solution of the vacuum Einstein equations, compatible with the initial data {'y,K). 

We note that in the original proof [29] of the well-posedness of the Einstein vacuum sys- 
tem, wave coordinates (then called "harmonic coordinates") are used. Wave coordinates 
correspond to a particular choice of the function Ta[g]- Generalized wave coordinates, 
which play an important role in numerical simulations of black hole collisions [52] corre- 
spond to other choices of J^a[g]- 

The techniques used to show that, for any given set of initial data (7, K) satisfying the 
constraint equations, there exists a unique maximal spacetime development of that data, 
are very different from those (just discussed) which are used to prove local existence. To 
prove that a unique maximal development exists, the idea is to consider the set A^[^ /^j of all 
spacetime developments of (7, ET), and then define a partial ordering ">" on A^[^^^], with 
(M^,^) > (M^,(7) if (M^,^) is an extension of (M^,^), up to diffeomorphism. One then 
shows that the ordered set {M.[^^k]^ >) has the properties needed to apply the Maximality 
Principle from set theory i. from which it follows that (TWr^ j^i, >) has a unique maximal 
element. One readily verifies that this maximal element is indeed a maximal spacetime 
development of the data (7, -RT). Details of this argument can be found in |13J and in 
Chapter 16 of [55]. 

There are two important ways in which one can generalize Theorem [TJ The first of 
these involves coupling in non-gravitational fields. It is not true that if a field theory has 



The Maximality Principle is equivalent to the Axiom of Choice. It is often referred to as "Zorn's 
Lemma, although some dispute the appropriateness of this label (see Chapter 16 of |55)). 
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a well-posed initial value problem in Minkowski spacetime, then the standard ("comma- 
to-semicolon"!^ coupling of that field theory to Einstein's theory also necessarily has a 
well-posed initial value problem. The standard Klein-Gordon vector field theory, with field 
equations D^DaW^ = mPW^, provides an example of a field theory which is well-posed in 
flat spacetime, yet appears to be ill-posed if coupled to Einstein's theoryri However, there 
are a number of field theories involving Einstein's theory coupled to non-gravitational fields 
that do have well-posed Cauchy problems. These include Einstein-Maxwell, Einstein- Yang- 
Mills, Einstein-DiracLj, certain forms of Einstein-scalar field theories, and certain forms of 
Einstein- fluid theories. For all of these, in addition to well-posedness, one can prove that 
for each set of initial data satisfying the appropriate set of constraints, there is a unique 
maximal development spacetime solution of the coupled system. 

The second way in which Theorem [1] can be generalized is by loosening the required 
degree of regularity. The original results, as well as the results we state here, presume 
that the initial data is smooth. Results of over forty years ago [28^ |35] show that well- 
posedness (as well as the existence of a unique maximal development) holds for initial 
data sets (7, K) with 7 contained in a local Sobolev spacq__| of index s > 5/2, and K in 
such a space with index s — 1. The work of Klainerman and Rodnianski [44J lowers the 
required regularity to s > 2, and the recent work of Klainerman, Rodnianski and Szeftel 
|45j indicates that (with certain restrictions) one can prove well-posedness for s = 2 as 
well. We note that the drive to achieve lower regularity well-posedness results is motivated 
both by the desire to understand weak solutions, and by the use of low regularity results 
as tools in understanding long-time behavior of solutions. 

Not stated as part of Theorem [T] above, but an important feature of Einstein's theory, 
is the fact that the quantities (7a6,-f^cd) induced on any Cauchy surface in a spacetime so- 
lution of Einstein's equations must satisfy the constraint equations. Hence, as one evolves 
a spacetime solution from a set of initial data, the evolution equations (|18p -(|19p effec- 
tively preserve the constraints. Since numerical implementation of the Cauchy problem 
inevitably introduces small errors, the constraints are not precisely preserved in a numeri- 
cally constructed solution. This appears to be a source of instability problems in numerical 
relativity. 



Roughly speaking, the "comma-to-semicolon" coupling of a given field theory to gravity involves re- 
placing partial derivatives with metric-compatible covariant derivatives in the Lagrangian for the given field 
theory, multiplying this Lagrangian by the metric volume density \/—detg, adding the result to the Einstein 
Lagrangian TZ^—detg, and then varying the summed Lagrangian with respect to the metric and the fields. 

It is generally difficult to prove that a system of partial differential equations does not have a well- 
posed Cauchy problem. However, for systems like the Klein-Gordon vector field theory with standard 
coupling to Einstein's theory, the presence of "derivative-coupling" terms in the equations cause standard 
analyses of hyperbolicity and well-posedness to be essentially unmanageable; one is led to strongly suspect 
ill-posedness. Some of these issues are discussed in [42j . 

The verification of well-posedness is straightforward for the Einstein-Dirac theory with commuting 
spinor fields. For Einstein-Dirac with anti-commuting spinor fields, and also for A'^ = 1 supergravity (which 
requires that the spinor fields anti-commute) , well-posedness holds, but in a subtle sense. See [4] for details. 

The Sobolev index s indicates (weak) L^ boundedness for order s derivatives of the indicated fields. 
See, e.g., [27] for definitions and properties of the Sobolev function spaces. 
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4. The conformal method of constructing and analyzing solutions of the 

constraint equations 

The crucial first step in building a spacetime solution of the Einstein equations via 
the initial value problem is to obtain a set of initial data which satisfies the constraints. 
Restricting our attention to the vacuum system for most of the discussion here, we presume 
that a fixed three-dimensional manifold E^ has been chosen, and we seek a Riemannian 
metric -^ab on S^ and a symmetric tensor Kd on S^ such that, together, {'^ab-,Kcd) satisfy 
the vacuum constraints (fT6l) - (fT71) . 

There are three somewhat different goals one might have in studying solutions of the 
constraint equations: i) the construction of physically interesting data sets, which may 
then be evolved into physically interesting spacetimes; ii) the parametrization (in terms 
of an appropriate function space) of the set of all solutions of the constraint equations; 
and iii) the systematic study of various mathematical and physical properties of solutions 
of the constraints, including (for asymptotically Euclidean data sets) conserved quantities 
such as global mass and global angular momentum. 

The most widely used analytical method for studying solutions of the constraints is the 
conformal method (together with the closely related conformal thin sandwich method). 
The conformal method is particularly adapted to the goal of finding a function space 
parametrization of the set of solutions of the constraints, the second goal listed above. It 
has also been used very effectively in the numerical construction of physically interesting 
initial data sets, especially those used to simulate black hole collisions (which are important 
for building gravitational signal templates). We discuss the conformal method here and the 
conformal thin sandwich method in Section [5l in Section [6l we discuss another approach 
for studying the constraints — gluing — which is especially effective as a tool for attaining 
the third goal listed above. 

As a PDE system for the twelve functions encompassed in ^ab{^) ^ind Kab{x), the four 
constraint equations (J16p -(jl7p constitute an underdetermined system. The idea of the 
conformal system is to split the initial data into two sets — the "free (conformal) data", 
and the "determined data"-in such a way that, for a specified choice of the free data, 
the constraint equations become a determined elliptical PDE system, to be solved for the 
determined data. There a number of ways to carry out this data split. We focus here 
on the "semi-decoupling split" (labeled in the early literature of the conformal method as 
"method A"; see [16)). which takes the following form: 

Free ("Conformal") Data: 

\ab ^ a Riemannian metric, specified up to conformal factor; 
(Tdfe - a divergence- free^^ (y^cTab = 0), trace- free (A°''o"ab = 0) 
symmetric tensor; 
r - a scalar field; 
Determined Data: 



In the conformal data, the divergence- free condition is defined using the Levi-Civita covariant derivative 
compatible with the conformal metric Aat- 
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(p - a positive definite scalar field; 
W"- - a vector field. 

For a given choice of the free data, the four equations to be solved for the four functions 
of the determined data take the form 

V^(Ll^)- = ^/^,T (25) 

Act> = lRcl)-ha"'^ + LW^^)ia„^n + LWmn)(p-' + ^T^(p\ (26) 

where the Laplacian A and the scalar curvature R are based on the Aaft-compatible covariant 
derivative Va, and where L is the corresponding conformal Killing operator, defined by 

iLW)ab := ^aWb + VbWa " ^A,feV„PF™. (27) 

Presuming that for the chosen conformal data one can indeed solve equations (j25p - (|26p for 
(j) and W, then the initial data set {'jab, Kcd) constructed via the formulas 

lab = (p'^Xab, (28) 

Kab = r\'Jab + LWab) + ^(p^XabT, (29) 

satisfies the Einstein constraint equations ()16p - ()17p . 

The four equations (j25p - (|26p . which we collectively refer to as the LCBY equations 
(since their derivation is based on the work of Lichnerowicz [56], Choquet-Bruhat, and 
York |16j). are readily obtained by substituting the formulas (|28p -(|29 p into the vacuum 
constraints (|16p-()17p. Two key identities play a major role in the derivation of (|25p - (|26p : 
The first is the formula for the scalar curvature of the conformally-transformed metric 
lab = (k'^^abi expressed in terms of the scalar curvature for Xab and derivatives of (/>: 

R{^) = (j)-*R{X) - 8Aa0. (30) 

We note that if we were to use a different power of (p as the conformal factor multiplying 
Aab, then this formula would involve squares of first derivatives of as well. The second key 
formula relates the divergences of a traceless symmetric tensor pab as calculated using two 
different covariant derivatives — V(^) (compatible with the metric 7) and Vm) (compatible 
with A): 

V^)P™b = r'V^)(0Vm6). (31) 

This formula dictates the choice of conformal scaling used in (|29p for the trace-free part of 
Kab] 3- different scaling would generally lead to the appearance of further V</) terms in the 
LCBY equations, which are to be avoided. 

Do the LCBY equations admit solutions ((/», W) for every choice of the conformal data 
(S^; Aab, (Ted, t)? It is easy to see that this is not the case: If we choose, for example, 
E^ = S^, Xab is the round metric on S^ , acd = everywhere, and r = 1 everywhere, then 
()25p takes the form 'Vm{LW)]^ = 0, which requires that LWab vanish everywhere. Equation 
(j26l) then takes the form A(/> = -gRcj) + j2^^ ■ Since the right hand side of this equation is 
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positive definite (recall the requirement that i;^ > 0), it follows from the maximum principle 
on closed (compact without boundary) manifolds that there is no solution. 

Since this example shows that solutions do not exist for every possible choice of the 
conformal data, we seek to determine exactly which sets of such data lead to a solution 
and which do not. This issue has been intensively studied for at least forty years, and while 
much is known, there remains much to be determined. Roughly speaking, for conformal 
data with constant mean curvature (CMC), we generally know for which sets of conformal 
data (S^; Aaf,,(Tcrf,r) solutions exist, and for which sets they do not. Uniqueness is well 
understood as well for CMC conformal data. For conformal data with nearly constant mean 
curvature, we know a number of classes of conformal data for which solutions exist, as well 
a number of classes for which they do not; there are, however, a number of unresolved 
cases. For conformal data with mean curvature far from constant, we know much less; 
however, there has been some progress recently in studying these sets of conformal data. 

To explain more specifically what is known and what is not known about solving the 
LCBY equations, it is very useful to classify conformal data sets into a wide collection of 
classes, based on a number of criteria, which include the following: 

• Manifold type and asymptotic conditions: Data on closed manifolds, asymptoti- 
cally Euclidean data, asymptotically hyperbolic data, asymptotically conical data, 
asymptotically cylindrical data, or data on manifolds with boundary. 

• Regularity conditions: Analytic data, smooth data, or data contained in specified 
Sobolev or Holder spaces. 

• Coupled non- gravitational fields: Vacuum Einstein, Einstein-Maxwell, Einstein- 
Dirac, Einstein-scalar, Einstein-fluid, or Einstein- Vlasov. 

It is beyond the scope of this article to discuss what is known for each of the several classes 
delineated by these criteria. Rather, we focus on what is known and what is not known 
for smooth data satisfying the vacuum Einstein constraints, either on closed manifolds or 
satisfying asymptotically Euclidean conditions. We include comments on some of the other 
classes at the end of this section. 

Within each class of conformal data, the key distinction is between those sets of data with 
constant mean curvature and those which have non-constant mean curvature. The mean 
curvature of any initial data set is given by the function trK := j'^^Kab = t (we use (|28l) - 
(]29p to calculate the last equality here), so we see that the CMC condition corresponds to 
choosing conformal data with constant r. Examining ()25p . we see that the CMC condition 
is important because constant r implies the vanishing of the right hand side of equation 
(j25p . which then resultCjin the vanishing of the tensor quantity {LW)ab- Consequently, for 
CMC conformal data the analysis of the LCBY equations reduces to that of the (decoupled) 
Lichnerowicz equation, which takes the form 

Ad) = ^Rcp - la^^amnr' + ^r^<l^\ (32) 



It follows from a straightforward analysis of the elliptic operator V • L that in essentially all cases, if 
m(I/W)r = 0, then LWab vanishes. 
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We note that for all of the non-gravitational fields coupled to gravity which are listed above, 
the first of the LCBY equations takes the form VmiLWy^ = ^(f>^Vg^T + Ja, where Ja 
depends on the non-gravitational conformal data I. and does not involve (/). Consequently, 
while LWab does not vanish for CMC data if non-gravitational fields are being considered, 
the LCBY equations do decouple, and as a result the determination if solutions to the 
LCBY equations exist depends essentially on the Lichnerowicz equation (with extra terms 
corresponding to the non-gravitational field^^^ 



4.1. CMC Data on Closed Manifolds. As noted above, for a set of constant mean 
curvature conformal data (S'^; Aafe,(Tcrf,r) for the vacuum Einstein equations, a solution to 
the LCBY equations exists (and a map to a solution of the constraint equations exists) if 
and only if there exists a solution to the Lichnerowicz equation ()32p . For this problem, we 
know exactly which sets of conformal data lead to solutions and which do not. To state 
these results (and to prove them), two features regarding conformal transformations of the 
conformal data are crucial. 

First, we note the Yamabe theorem [56] for Riemannian metrics on closed manifolds, 
which states that every such metric can be conformally transformed to a metric of constant 
scalar curvature; further, for each metric the sign of that constant scalar curvature is 
unique. Hence, the set of Riemannian metrics on a given manifold S^ are partitioned into 
three "Yamabe classes" y~^,y^, and y~ , depending on that sign. 

Second, we readily verify that the Lichnerowicz equation is conformally covariant in the 
following sense: If is a solution to the Lichnerowicz equation for a set of conformal data 
(E^; Aah, o'cd, t), then for any smooth positive function 0, the function cp = 9~^ (pis a solution 
to the Lichnerowicz equation for the conformal data (S^; 0'^Xab, ^"^^cd; ''")• Combining this 
result with the Yamabe theorem, we see that a solution to the Lichnerowicz equation exists 
for (T,^; Xab,o'cd,T) if o-nd only if there exists a solution for the conformally transformed 
data iTi^]6^Xab-,0~'^acdiT) with R[9'^Xab\ = +1,0 or —1, depending upon the Yamabe class 
of Xab- 

Besides the Yamabe class of Xab, two aspects of a given set of conformal data determine 
whether or not a solution to the Lichnerowicz equation exists. One of them is whether or 
not the constant r is zero or not. In the table of results below, we label these alternatives 
r = or r 7^ 0. The other is whether or not the function acd<^'^'^ is identically zero (on S^) 
or not. We label these alternatives o" = or a ^ 0. Using "Y" to indicate that solutions 
exist for conformal data sets in a certain class and using "N" to indicate that they do not, 
we summarize the results in the following table: 



This effect occurs as a consequence of the choice one makes for the conformal rescahng of the non- 
gravitational fields. Alternative choices could be made which would lead to Ja depending on </!>, but there 
is little motivation for making such choices. 

For example, in the case of the Einstein-Maxwell theory, the Lichnerowicz equation picks up a term 
of the form —p(f>~^, where p is quadratic in the conformal electric and magnetic fields. 
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We note that for all those data sets such that solutions exist, except for that class of data 
with X ^ y^, with (7 = and with r = 0, the solutions are unique. In the latter (somewhat 
trivial) case, any (positive) constant cp is a solution. 

The results summarized in the table here were to a large extent proven in the 1970's by 
Yvonne Choquet-Bruhat, James York, and Niall O'Murchadha; see [16] for a discussion of 
this early work. The complete proof, including two cases not handled by the previous work, 
appears in [36] . There it is shown that the "No" cases can all be proven using the Maximum 
Principle, stated in the following form: On a closed manifold, the equation A(j) = f[x,(j)), 
with /(x, (p) either non-vanishing and non-positive or non-vanishing and non-negative, has 
no solution. It is also shown there that the "Yes" cases can all be proven using the Suh and 
Super Solution Theorem, which can be stated as follows: If there exist a pair of positive 
functions (/>_(. > (/>_ which satisfy the inequalities 



A0_ < -R(t) 



1 



mn 



+ 



1 

12' 



^2 i5 



and 



1 1 

> -R(p+ - -a'' 



(Jr, 



' +12^^ 



(33) 



(34) 



then there exists a solution (j) of the Lichnerowicz equation ([32]), with (/)+></>></>-. For 
some of the six "Yes" cases (specifically those with the metric in the negative Yamabe 
class), constant sub and super solutions are easily found. For others, non-constant sub and 
super solutions are needed, and are not so easily found. 

We can use the results summarized in the table, along with the conformal covariance 
stated above and a certain scaling invariance, to parametrize the set of CMC solutions of 
the vacuum constraints on a chosen closed manifold E^. We note first that every solu- 
tion {'jabjKab) of the constraint equations (fT6]) - (fT7l) can be obtained using the conformal 
method, since one can choose Xab = Jab, t = K^, and acd equal to the divergence-free trace- 
free projection i of Kcd', with this conformal data, W^ = and <p = 1 satisfy the LCBY 
equations, leading us back to (7afe, Kab)- Next, we note that the conformal covariance result 
implies that for any positive function 9 the two sets of conformal data (S^; Xab, (^cd, t) and 
{Tj^;9^Xab,0~'^acd,T) lead to related solutions of the Lichnerowicz equation and thence to 
identical solutions {'jab, Kcd) of the constraint equations ([T6l) - (fT7l) . The scaling result [5] 
states that for a positive constant A, the two sets of CMC conformal data (S^; Aafe,o"cd,T) 
and (S^; A^Xab, Aacd, A~^t) admit the same solution 4> of the Lichnerowicz equation. These 
lead to different solutions {'yab,Kcd) and {A^jab, AKcd) of the constraint equations; but 



As discussed in 16], such a projection always exists 
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these solutions are related by an essentially trivial scaling. As a consequence of these con- 
siderations, we determine that for any chosen closed manifold S^, the set of all conformal 
data sets marked "Yes" in the table, quotiented out by the conformal covariance equiva- 
lence and by the scaling equivalence, provides a faithful (bijective) parametrization of the 
space of CMC solutions of the vacuum constraints on S^. 

These results for CMC solutions of the vacuum constraints are essentially replicated 
if one considers the Einstein-Maxwell, Einstein- Yang-Mills, Einstein-Dirac, or Einstein- 
fiuid 1 constraint equations [38]. On the other hand, the results for the Einstein-scalar field 
theories with certain types of field theory potentials are not nearly as complete. Details of 
the difficulties that can arise are discussed in |14j. 

4.2. Asymptotically Euclidean CMC Data. The decoupling of the LCBY equations 
which results from the assumption of constant mean curvature does not depend upon the 
topology of S^ or the asymptotic conditions of the conformal data. Hence, as for the closed 
manifold case discussed above, the analysis of the conformal method for asymptotically Eu- 
clidean data (as well as for asymptotically hyperbolic data or any other chosen asymptotic 
condition) reduces to the study of the solvability of the Lichnerowicz equation. We note 
that the fall-off conditions which define a set of data to be asymptotically Euclidean (see, 
for example, [15j for a precise statement of these fall-off conditions) require that the mean 
curvature, if constant, be zero. Hence asymptotically Euclidean data which is CMC must 
be maximal; i.e., K^ = r = 0. 

The criterion for the Lichnerowicz equation to admit a solution for asymptotically Eu- 
clidean conformal data is independent of acd, and depends only the metric. As proven 
by Cantor [12], the Lichnerowicz equation admits a solution for a given set of asymptoti- 
cally Euclidean conformal data if and only if the metric admits a conformal transformation 
which results in the scalar curvature vanishing everywhere. Such metrics are called "Yam- 
abe positive", and Brill and Cantor [TT] have also shown (with a correction pointed out 
by Maxwell [17]) that an asymptotically Euclidean metric is Yamabe positive in the sense 
just described if and only if for every non-vanishing compactly supported function / on 
S^, the inequality 

,nf J"'l^^l'„t,f)^^>0. ,35) 

holds. 

We note that the same criterion for the solvability of the Lichnerowicz equation holds 
for asymptotically Euclidean data for the Einstein-Maxwell, Einstein- Yang-Mills, Einstein- 
Dirac, Einstein-fiuid, and other such theories involving non-gravitational fields coupled to 
Einstein's theory.. 



For the Einstein-fluid constraints, there are consistent choices of conformal scaling of the fluid flelds 
which lead to a more difficult set of LCBY equations, analogous to those arising for certain Einstein-scalar 
theories. Conformal scalings can, however, always be chosen in a way which avoids these problems [38j . 
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4.3. Near-CMC Data. If we drop the non-CMC condition on the choice of the conformal 
data, then Vr does not vanish, and we must deal with the fully coupled LCBY system. 
While not much is known about the solvability of the LCBY system for general sets of non- 
CMC conformal data, for those sets of data with | Vr| small in some appropriate sense, one 
can in many cases determine whether or not solutions exist. 

One way to analyze the coupled system is to use the iterated Gummel method, which 
replaces ([25|) - ([26j) by the sequence of semi-decoupled PDE systems 






S'^Cn-l) 



1 /„cd 



+ LW[^)){acd + LWin)cd) 
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(36) 
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The idea is to (i) choose an initializing value for (J)q] (ii) show that a sequence (0(n);W^(n)) 
of solutions to (i36]l - (f37|l exists; (iii) prove that there are uniform upper and lower bounds 
for the sequence ((/<(„), Ty(„)); and (iv) use those upper and lower bounds together with a 
contraction mapping argument to show that the sequence {(puA, VF(n)) converges uniformly 
to a limit ((/), W), which solves the LCBY equations. 

The assumption that |Vt| is small plays a crucial role in carrying out steps (iii) and (iv). 
It does this because, as a consequence of the elliptic equation ([25]) and its sequential analog 
(j36p . the norm of LW is controlled by | Vt|0^. Thence the term involving the square of LW 



in (1261). which contains factors of 



tfi\2 



times 



-I.e., 



—competes with the t (jr term. 
To ensure that t'^<J)^ (which has the favorable sign) wins this competition and to thereby 
retain the control of 4> that one has in the CMC case, it is sufficient to require that Vrr be 
sufficiently small 

Analyses of the sort discussed above have been carried out in |40] [3]. Combining these 
results with others proven in [37] and [25], one has the following table which lists, for each 
of twelve classes of near-CMC conformal data on closed manifolds, whether solutions to 
the LCBY equations are known to exislo ("Y"), are known to not exist ("N"), or are not 
known either way ("?"). In this table, the Yamabe classes and their labels {3^~'',3^'^,3^~} 
are the same as used above for the CMC table, and the labels o" = and o" ^ denoting 
whether or not \(t\ is identically zero or not are also the same as for the CMC table. For 
the function r, the notation t'^ > indicates those sets of conformal data with r (presumed 
smooth) having no zeros, while r^ ^ labels those sets of conformal data with r allowed 
to have zeroes. 
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Uniqueness holds as well in all of these "Y" cases. 



THE INITIAL VALUE PROBLEM IN GENERAL RELATIVIT^ 17 

Not surprisingly, the classes which cause the most difficulty are those in which r has 
zeroes; we have a fairly complete understanding of whether or not solutions exist for the 
classes in which r is bounded away from zero. 

Near-CMC conformal data sets which are asymptotically Euclidean have also been stud- 
ied. Techniques similar to those discussed briefly here show [?] that the criterion for the 
existence of solutions to the LCBY equations for near-CMC asymptotically Euclidean con- 
formal data sets is very similar to that for CMC data. 

4.4. Far-CMC Data. Without either the CMC or the near-CMC condition imposed, the 
analysis of the LCBY equations is considerably more difficult. For most sets of conformal 
data satisfying neither of these conditions (we use the label "far-CMC" both for initial 
data sets (S^; 7^6, Kcd) and for conformal data sets (S^; Aafe,o"cd,r) which are not CMC 
and satisfy no smallness condition on K^ or on r), it is not known whether or not solutions 
exist. 

There has, however, been recent work which begins to explore whether or not solutions 
exist for far-CMC data, and also explores their multiplicity. We briefly discuss three of 
these recent works here. 

The work of Hoist, Nagy, and Tsogtgerel [Mj together with the important follow-up 
work of Maxwell [58] effectively swaps the near-CMC condition of smallness of |Vr| for a 
condition requiring the smallness of \(Jcd\- More specifically, their combined work shows that 
the LCBY equations admit a solution for a given set of conformal data on a closed manifold 
if the following conditions hold: (i) Xab £ y^', (ii) ^ab does not admit a conformal Killing 
field; (iii) \crcd\ is sufficiently small; and (iv) acd is not identically zero. We emphasize the 
fact that these conditions include no restriction on the mean curvature function r. We also 
note that while these conditions guarantee existence, they tell us nothing about uniqueness. 

Besides showing that a small but interesting class of far-CMC conformal data is mapped 
to solutions of the constraint equations, this work introduces a new analytical tool to the 
study of the conformal method. The proof that the LCBY equations admit solutions for 
sets of conformal data of this sort described above relies on expressing solutions of the 
LCBY equation formally as fixed points of a map 

^ : <A ^ A-i [ii?0 - i(^ + L{(V • L)-\^<p'VT)}fr' + ^r'] , (38) 

and proving that such fixed points exist using Schauder compactness. 

Schauder techniques also play a role in the work Dahl, Gicquaud and Humbert [?]. 
They show that for certain classes of conformal data — data on a closed manifold S^ with 
the metric (any Yamabe class) admitting no conformal Killing fields, with r bounded away 
from zero, and with |(Tp not identically zero if Xab G y~ — solutions of the LCBY equations 
exist if the equation 

Ay" = |Ly|-W (39) 

r 

does not admit a solution. Thus, from this perspective, one proves the existence of LCBY 

solutions for a given set of conformal data by showing that (I39p (called the limit equation 

by the authors of ^25j) admits no solution. This approach does not directly prove the 
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existence of solutions of the LCBY equation for sets of far-CMC conformal data, but it 
can in principle be used to do this. In [31j, it shown that an analogous limit equation-type 
result holds for asymptotically hyperbolic conformal data sets. 

It is interesting that the derivation of these limit equation results relies on the study 
of solutions of a family of "subcritical LCBY^ equations" in which equation ()26p is left 
unchanged, but ([25|) is replaced by 

Vm{LW)^ = ^</.(^-^)v,r (40) 

It is much easier to prove that solutions to the system (I40p - (l26p . with e > 0, exist, than 
to prove the existence of solutions to the LCBY equations. The idea then is to study the 
e — 7- limit of solutions of the LCBY^ system, and use the analysis of these limits to shed 
light on the LCBY system. Such an analysis leads to the limit equation results. 

The third work concerning solutions of the LCBY equations for far-CMC conformal data 
which we discuss here is that of Maxwell [49] , in which he studies a very simple class of 
planar symmetric conformal data and finds somewhat surprising results. The conformal 
data sets he considers are characterized as follows: (i) T?' = T^; (ii) the metric is flat; 
(iii) acd has planar (T^) symmetry, and therefore, as a consequence of the trace- free and 
divergence- free conditions, is a matrix with constant entries (two free constants); and (iv) 
r is a step function, with a pair of discontinuities. Working with these simple data sets, 
which are parametrized by a small set of constants. Maxwell can directly (numerically) 
construct solutions, if they exist. He finds that for certain ranges values of the parameters, 
no solutions exist; for other ranges, multiple solutions exist; and finally for a third range, 
unique solutions can be obtained. 

It is very intriguing to consider whether these results of Maxwell generalize to conformal 
data sets without any discontinuities, and to data sets with less imposed symmetry. In 
any case, these results suggest that the behavior of the conformal method for far-CMC 
conformal data sets could be interesting and complicated. 

5. The Conformal Thin Sandwich Method as an Alternative to the 

Conformal Method 

The conformal method has proven to be a remarkably useful tool for generating and 
parametrizing and analyzing solutions of the Einstein constraint equations. It does, how- 
ever, have some minor drawbacks: a) The conformal data is somewhat remote from the 
physical data, since the conformal factor changes the physical scale on different regions of 
space, b) While casting the constraints into a determined PDE form has the advantage of 
producing PDEs of a relatively familiar (elliptic) form, one does give up certain flexibilities 
which are inherent in an underdetermined set of PDEs. c) In choosing a set of conformal 
data, one has to first project out a divergence-free trace-free tensor field {(Jed)- d) While 
the LCBY system is conformally covariant in the sense discussed above in Section [4. II for 
CMC conformal data, this is not the case for non-CMC conformal data. 

The last two of these problems can be removed by modifying the conformal method in 
a way which York ^57j has called the "conformal thin sandwich" (CTS) approach. The 
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basic idea of the conforinal thin sandwich approach is essentially the same as that of the 
conformal method. There are, however, two important differences. First, the CTS free data 
sets are larger than the free (conformal) data sets of the conformal method in the following 
sense: Like the conformal data sets, the CTS data sets include a conformal metric Xab and 
a mean curvature scalar r. In addition, the CTS data sets include a trace-free tensor Ucd 
to replace the divergence-free trace-free tensor acd of the conformal data sets, plus an extra 
scalar field rj. Second, after solving the following set of CTS equations (analogous to the 
LCBY equations) 

V„((2r?)-i(LX))- = h<^V^r + V^{{2r^)-'UT), (41) 

A$ = ii^cI>_i([/™- + Ly-")(C/„„ + Ly„„)$-^ + lr2$^ (42) 

for the conformal factor <1> and the vector field V^, one constructs not just the initial data 
(7afe, Kcd), but the lapse A^ and the shift M" as well: 

lab = ^^Xab (43) 

Kab = <^-\-Uab + LYab) + ^^^XabT (44) 

N = ^^7] (45) 

M" = Y". (46) 

It is clear that in using the CTS approach, one need not project out a divergence-free 
part of a symmetric trace- free tensor. As well, one also readily checks that the CTS method 
is conformally covariant in the sense discussed above: the initial data (7a6,-f^cd) and the 
lapse and shift {N, M") generated from the CTS data set {Xab, Uab, t, rj) and from the CTS 
data set {6'^Xab,0~'^Uab,T,9^ri) are identical. Furthermore, since the mathematical form 
of equations (j4ip -()42p is very similar to that of (|25p - (j26p . the solvability results for the 
conformal method can be essentially carried over to the CTS approach. 

There is, however, one problematic feature of the conformal thin sandwich approach. 
The problem arises if we seek CMC initial data with the lapse function chosen so that 
the evolving data continues to have constant mean curvaturq_j. In the case of the con- 
formal method, after solving (I25p - ()26p to obtain initial data (7ab,-f^cd) which satisfies the 
constraints, one achieves this by proceeding to solve a linear homogeneous elliptic PDF for 
the lapse function. One easily verifies that solutions to this extra equation always exist. 
By contrast, in the CTS approach, the extra equation takes the form 

A($^r/) = -'^'^rjR + -($7?)^i(C/ - LY)"^ + ^^y^V^r - $^ (47) 

8 2 



which is coupled to the system ()4ip -()42 p . The coupling is fairly intricate; hence little is 
known about the existence of solutions to the system, and it has been seen that there 
are problems with uniqueness. These difficulties do not arise, of course, if one makes no 
attempt to preserve the constant mean curvature condition. 



on 

Such a gauge choice is often used in numerical relativity. 
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6. Gluing Solutions of the Constraint Equations 

Both the conformal method and the conformal thin sandwich method are procedures for 
generating initial data sets which satisfy the Einstein constraint equations from scratch. 
The gluing procedures, which we discuss here, produce new solutions of the constraint 
equations by combining existing solutions. While the gluing procedures have not yet turned 
out to be as useful as the conformal method and CTS method for the practical generation 
of physical interesting initial data sets, they have proven to be very effective for certain 
applications and for settling certain conjectures. We outline some of these applications 
below, after describing the two gluing procedures which have been developed, and how 
they work. 

The asymptotic exterior gluing, developed by Corvino and Schoen [231 ^^ works as 
follows. We presume that {T,'^,"fab,Kcd) is an asymptotically Euclidean initial data set 
which satisfies the Einstein constraints, and also satisfies certain asymptotic conditions (as 
specified in |24]). For any compact region H^ C S^ for which S^ \ H^ = M'^ \ B^ (where 
B^ is a ball in M^), there is a smooth asymptotically Euclidean solution of the constraints 
on T,^ which is identical to the original solution on H^ C T,^, and is identical to Cauchy 
data for the Kerr solution on S'^ \ S'^ for some H^ C S^. In words, this technique allows 
one to smoothly glue any interior region of an asymptotically Euclidean solution to an 
exterior region of a slice of a Kerr solution. There is generally a transition zone between 
the interior chosen region and the exterior Kerr region S^ \ H which is unknown, but is a 
solution of the constraint equations. We note that for asymptotically Euclidean solutions 
of the constraints with K^d = 0, this method glues any interior region to an exterior region 
of a slice of the Schwarzschild solution. 

In some sense, the Corvino-Schoen asymptotic exterior gluing result is very surprising. 
// the constraint equations were a determined elliptic system, one would not expect to 
be able to smoothly glue two solutions together like this, even with a transition region 
(satisfying the constraints). The key to proving that asymptotic exterior gluing indeed 
works is the exploitation of the underdetermined character of the constraints as a PDE 
system. Some of the ideas developed in the Corvino-Schoen work have also proven to be 
useful for localizing metric deformations as solutions of the constraints, as shown in the 
work of Chrusciel and Delay |21j . 

The other gluing procedure, connected sum gluing, was developed first by Isenberg, 
Mazzeo and Pollack [391 with further work done together with Chrusciel p2| and with 
Maxwell [38] . The idea here is to start with a pair of solutions of the (vacuum) constraints 
{T,f,ji,Ki) and (S2,72,^2) and to choose of a pair of points pi G Sf and p2 G ^l, one 
point contained in each solution. Based on these solutions, connected-sum gluing produces 
a new set of initial data (Sq_2n,7(i_2))-^(i-2)) with the following properties: i) 5](i_2) is 
diffeomorphic to the connected sunoSf^Sl; ii) (S?^ ,,-,,7(i_2)) -?^(i-2)) is a solution of the 



The connected sum of these two manifolds is constructed as follows: First we remove a ball from each 
of the manifolds Ef and E2. We then use a cylindrical bridge S^ x I (where / is an interval in R^) to 



connect the resulting S^ boundaries on each manifold. 
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constraints everywhere on Sq_2\; iii) On that portion of 5]q_2\ which corresponds to Sf \ 
{ball around pi}, the data (7(i_2))-f^(i-2)) is isomorphic to {'yi,Ki), with a corresponding 
property holding on that portion of S^ which corresponds to E^ \ {ball around p2}- 

Connected sum gluing can be carried out for fairly general sets of initial data. The 
sets may be asymptotically Euclidean, asymptotically hyperbolic, specified on a closed 
manifold, or indeed anything else. The only condition that the data sets must satisfy is that, 
in sufficiently small neighborhoods of each of the points at which the gluing is to be done, 
there do not exist nontrivial solutions ^ to the equation DQ*, ^^^C = Oj where DQ*, ^n is 

the operator obtained by taking the adjoint of the linearized constraint operatoCj. In work 
by Beig, Chrusciel and Schoen [9], it is shown that this condition (sometimes referred to 
as "No KIDs" , meaning "no (localized) Killing initial data" ) is indeed generically satisfied. 

The details of the proof that connected sum gluing can be carried out as generally as 
described above are beyond the scope of this paper; see [22] along with the references 
cited in that work for a complete discussion. We do wish to note three features of the 
proof: First, the proof is constructive in the sense that it outlines a systematic, step-by- 
step mathematical procedure for doing the gluing. In principle, one should be able to 
carry out the gluing procedure numerically. Second, connected sum gluing relies primarily 
on the conformal method, but it also requires the use of a non-conformal deformation 
(dependent on the techniques of Corvino and Schoen, and of Chrusciel and Delay), so as 
to guarantee that the glued data is not just very close to the given data on regions away 
from the connecting bridge, but is indeed identical to it. Third, while the Corvino-Schoen 
asymptotic exterior gluing has not yet been proven to work for solutions of the constraints 
with source fields, connected sum gluing (up to the last step, which relies on Corvino- 
Schoen) has been shown to work for most matter source fields of interest [38] . It has also 
been shown to work for general dimensions greater than or equal to three. 

As noted above, while gluing has not seen wide-spread use as a procedure for producing 
physically interesting initial data sets, it has proven to be very valuable for a number of 
applications. We note a collection of these applications here: 

(1) Spacetimes with Regular Asymptotic Structure: Until recently, it was not known 
whether there is a large class of spacetime solutions of the Einstein equations which 
admit the conformal compactification and consequent asymptotically simple struc- 
ture at null and spacelike infinity characteristic of the Minkowski and Schwarzschild 
spacetimes. Using asymptotic exterior gluing, together with Friedrich's analyses of 
spacetime asymptotic structures and arguments of Chrusciel and Delay |20|, one 
can produce such a class of solutions. 

(2) Initial Data for the Gravitational N-Body Problem: To model the physics of a sys- 
tem consisting of N chosen astrophysical bodies interacting gravitationally, it is 
important to be able to construct initial data sets which solve the Einstein con- 
straints and which accurately model the bodies of interest, their initial placement. 



If a solution to this equation does exist on some region A £ E^, it follows from the work of Moncrief 
that the spacetime development of the data on A admits a nontrivial isometry. 
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and their initial momenta, all in a single asymptotically Euclidean space. Chru- 
sciel, Corvino, and Isenberg [18] have used gluing techniques to show that for any 
chosen set of N asymptotically Euclidean solutions of the constraints representing 
black holes, stars, or other astrophysical objects of interest, one can construct a 
new asymptotically Euclidean solution which includes interior regions of these N 
chosen solutions, placed as desired (so long as the distances between the bodies are 
sufficiently large) and with the desired relative momenta. 

(3) Adding a Black Hole to a Cosmological Spacetime: Although there is no clear 
established definition for a black hole in a spatially compact solution of Einstein's 
equations, one can glue an asymptotically Euclidean solution of the constraints to 
a solution on a compact manifold, in such a way that there is an apparent horizon 
on the connecting bridge. Studying the nature of these solutions of the constraints, 
and their evolution, could be useful in trying to understand what one might mean 
by a black hole in a cosmological spacetime. 

(4) Adding a Wornihole to Your Spacetime: While we have discussed connected sum 
gluing as a procedure which builds solutions of the constraints with a bridge con- 
necting two points on different manifolds, it can also be used to build a solution 
with a bridge connecting a pair of points on the same manifold. This allows one 
to do the following: If one has a globally hyperbolic spacetime solution of Ein- 
stein's equations, one can choose a Cauchy surface for that solution, choose a pair 
of points on that Cauchy surface, and glue the solution to itself via a bridge from 
one of these points to the other. If one now evolves this glued-together initial data 
into a spacetime, it will likely become singular very quickly because of the collapse 
of the bridge. Until the singularity develops, however, the solution is essentially as 
it was before the gluing, with the addition of an effective wormhole. Hence, this 
procedure can be used to glue a wormhole onto a generic spacetime solution. 

(5) Removing Topological Obstructions for Constraint Solutions: We know that every 
closed three dimensional manifold S'^ admits a solution of the vacuum constraint 
equations. To show this, we use the fact that S^ always admits a metric F of con- 
stant negative scalar curvature. One easily verifies that the data {'j = T,K = T) 
is a CMC solution. Combining this result with connected sum gluing, one can 
show that for every closed S^, the manifold S^ \ {p} admits both an asymptoti- 
cally Euclidean and an asymptotically hyperbolic solution of the vacuum constraint 
equations. 

(6) Proving the Existence of Vacuum Solutions on Closed Manifolds with No CMC 
Cauchy Surface: Based on the work of Bartnik [6j one can show that if one has 
a set of initial data on the manifold T^^T^ with the metric components even- 
reflective across a central sphere and the components of K odd-reflective across 
that same central sphere, then the spacetime development of that data does not 
admit a CMC Cauchy surface. Using connected sum gluing, one can show that 
indeed initial data sets of this sort exist 1221. 
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7. Comments on the Long-Time Evolution of Spacetime Solutions of 

Einstein's Equations 

Once an initial data set satisfying the Einstein constraint equations has been obtained, 
one can evolve it into a spacetime satisfying the Einstein field equations. As guaranteed by 
the work discussed in Section^ there is a unique globally hyperbolic spacetime development 
of this initial data which contains (up to diffeomorphism) any other developments of the 
same set of data. 

What do we know about the long-time properties of these maximal developments? The 
Hawking-Penrose singularity theorems |33J tell us that (among spacetimes with a compact 
Cauchy surface), in one or the other direction in time, such developments "generically" 
become causally geodesically incomplete, which means that there are causal geodesies in 
the spacetime which do not extend to infinite affine parameter length. This property of 
causal geodesic incompleteness is, however, consistent with a wide variety of spacetime 
behavior, including curvature blowup, Cauchy horizon formation, and various topological 
anomalies [l]. 

One of the intriguing questions concerning spacetime developments is which of these 
behaviors — curvature blowup, Cauchy horizon formation, or something else — is expected 
to occur generically among those spacetimes which are causally geodesically incomplete. 
Penrose J5T] has conjectured that curvature blowup is the generic behavior. This conjecture 
has been labeled the Strong Cosmic Censorship conjecture (SCCjfj and it is viewed by 
many as one of the central questions concerning the evolutionary behavior of solutions of 
Einstein's equations. 

It is well-known that there are infinite dimensional families of solutions [50j which have 
bounded curvature and develop Cauchy horizons. The existence of these solutions does not 
disprove SCC. To formulate and study the SCC conjecture carefully, one needs to define 
the notion of "generic" in terms of the topology of the space of constraint-satisfying initial 
data sets on a fixed three dimensional manifold S^, and then determine which sets of initial 
data evolve into a spacetime with unbounded curvature, and which do not. 

Strong Cosmic Censorship, formulated this way, has been proven for certain families of 
solutions, most notably (by Ringstrom ^53j) for the Gowdy familjo which is characterized 
by the existence of a T^ isometry group, and vanishing "twists". Numerical and other 
formal evidence strongly suggest that it is true for a wider class of spacetimes, with smaller 
isometry group f4T] . Proving or disproving this remains a major challenge. 

The Strong Cosmic Censorship conjecture concerns generic behavior among those space- 
times which are causally geodesically incomplete. Distinct from this issue, and also of very 
significant interest is the question of which initial data sets evolve into spacetimes which 



The Strong Cosmic Censorships conjecture does not imply, and is not implied by, the Weak Cosmic 
Censorship conjecture, which concerns the generic formation of an event horizon around a singularity which 
forms in an asymptotically flat solution of Einstein's equations. 

The Gowdy family of solutions is introduced and characterized in 32 and is studied extensively in 
[19j . SCC is proved for the Gowdy spacetimes in [53] . 
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extend an infinite (proper) time into the future and/or the past, and which do not. A com- 
plete answer to this question appears to be beyond our current mathematical capabilities. 
However, as a small but very significant step towards answering this question, a number 
of researchers have focussed on the issue of the stability — in terms of long time existence 
and structure — of solutions which exist for infinite proper time. 

The landmark stability result in general relativity is the proof by Christodoulou and 
Klainerman p.7j of the stability of Minkowski spacetime. They show, using energies based 
on the Bel-Robinson tensor to measure initial data perturbations, that the spacetime de- 
velopments of initial data sets which are sufficiently close to Minkowski initial data do 
extend an infinite proper time into the future and into the past. Moreover, they show that 
these developments have the same asymptotic spacetime structure as Minkowski space- 
time. These results have been extended to allow electromagnetic as well as gravitational 
initial data perturbations [10] and have also been strengthened in terms of the nature of 
the asymptotic structure which is shown to be stable [TO]. 

If Minkowski spacetime is stable, one might logically proceed to consider if this is also 
the case for Schwarzschild spacetimes. However, since one knows that a small perturbation 
of Schwarzschild initial data which adds angular momentum will evolve into a Kerr solution 
rather than a Schwarzschild solution, one is led to consider the stability of Kerr spacetimes 
instead. The determination of whether or not Kerr spacetimes are stable is currently one 
of the most active areas of research in mathematical relativity. A recent report on the 
research directed towards this goal appears in |26j . 

We note that the stability of other spacetimes has been established: Friedrichs pO] has 
shown that DeSitter spacetime is stable, Andersson and Moncrief [2] have shown that Milne 
spacetime is stable, and Ringstrom [53] has shown that certain solutions of the Einstein- 
scalar field equations with accelerating expansion are stable. 
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